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ABSTRACT 

We present an absorption line analysis of the Lyman limit system (LLS) at z « 3.55 in our Mag- 
ellan/MIKE spectrum of PKS2000— 330. Our analysis of the Lyman limit and full H I Lyman series 



constrains the total H I column density of the LLS (iVni = 10 



18.0±0.25 



cm for &hi > 20 km s ) and 



also the Nm values of the velocity subsystems comprising the absorber. We measure ionic column 
densities for metal-line transitions associated with the subsystems and use these values to constrain 
the ionization state (> 90% ionized) and relative abundances of the gas. We find an order of magni- 
tude dispersion in the mctallicitics of the subsystems, marking the first detailed analysis of metallicity 
variations in an optically thick absorber. The results indicate that metals are not well mixed within 
the gas surrounding high z galaxies. Assuming a single-phase photoionization model, we also derive 
an Vn-weighted metallicity, < [Si/H] >= — 1.66 ± 0.25, which matches the mean metallicity in the 
neutral ISM in high z damped Lya systems (DLAs). Because the line density of LLSs is over 10 x 
higher than the DLAs, we propose that the former dominate the metal mass-density at z ~ 3 and 
that these metals reside in the galaxy/IGM interface. Considerations of a multi-phase model do not 
qualitatively change these conclusions. Finally, we comment on an anomalously large O /Si + ratio 
in the LLS that suggests an ionizing radiation field dominated by soft UV sources (e.g. a starburst 
galaxy). Additional abundance analysis is performed on the super-LLS systems at z w 3.19. 
Subject headings: large-scale structure of universe — quasars: absorption lines — intergalactic medium 



1. INTRODUCTION 



Lyman limit systems (LLSs) are the 'clouds' along 
quasar sightlines that have sufficient H I column den- 
sity Nm to be optically thick (r > 1) at the Lyman 



limit, i.e. 7Vhi > 



This definition sepa- 



rates LLSs from the Lya forest lines (clouds with typical 
Nm < 10 15 cm~ 2 ) which trace the intergalactic medium 
(IGM; Rauch 1998) and the damped Lya systems (DLAs; 
Nm > 2 x 10 20 cm -2 ) that trace high redshift proto- 
galaxies (Wolfe et al. 2005). On the basis of their inter- 
mediate H I surface density, one may expect the LLSs 
to represent the physical interface between the IGM and 
high z galaxies. Studies of the LLSs, therefore, may con- 
strain the accretion history of galaxies and the transport 
of material back to the IGM. 

A key tracer of these processes is the metals inter- 
mixed with the gas. These are presumably formed in 
the deep potential wells of galaxies, imported to the in- 
terstellar medium by stellar winds and supernova explo- 
sions, and (possibly) carried to greater distances by tidal 
disruption and/or star- formation and AGN feedback pro- 
cesses. With the advent of the 10m Keck I telescope cou- 
pled with the HIRES echelle spectrometer (Vogt et al. 
1994), observers demonstrated that a significant frac- 
tion of the IGM is enriched in heavy elements of C, Si, 
and O (Tytler et al. 1995; Cowie et al. 1995). More re- 
cent studies have established the variation of metallicity 
with density and redshift (Schaye et al. 2003; Simcoe 
et al. 2004; Aguirre et al. 2008) and demonstrate a me- 
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dian metallicity of roughly 1/1000 solar. These obser- 
vations constrain the integrated enrichment by galaxies 
of their surrounding medium (e.g. Aguirre et al. 2001). 
Studies of the metallicity distribution and evolution in 
the damped Lya systems (Pettini et al. 1994; Prochaska 
et al. 2003), meanwhile, constrain the star- formation his- 
tory of high z galaxies and the processes of ISM enrich- 
ment (Dessauges-Zavadsky et al. 2004; Johansson & Ef- 
stathiou 2006). These studies demonstrate that the ISM 
of most high z galaxies is metal-poor (w 1/30 solar), that 
none have a metallicity less than 1/1000 solar, and that 
a wide dispersion (lOOx) exists from galaxy to galaxy 
(Prochaska et al. 2003). 

Although precise measurements for over 100 damped 
Lya galaxies have been acquired (Prochaska et al. 2003; 
Ledoux et al. 2006; Prochaska et al. 2007), constraints 
on metallicity variations within single galaxies are very 
limited. The metal-line transitions generally show ab- 
sorption by multiple components spanning « 100 km s _1 
(Prochaska & Wolfe 1997), presumably by gas dis- 
tributed on galactic (kpc) scales. The Lyman series, 
however, is in general too saturated to yield the H I 
column densities of the velocity components identified 
in the metal-line transitions. Therefore, these data pro- 
vide little constraint on metallicity gradients along the 
sightline even as a function of velocity much less posi- 
tion. The studies to date have focused instead on rel- 
ative abundance ratios within DLAs which best trace 
variations in the gas density, depletion, and/or ioniza- 
tion state (Prochaska & Wolfe 1996; Lopez et al. 2002; 
Prochaska 2003; Dessauges-Zavadsky et al. 2006). 

In contrast to DLAs, the higher order Lyman series 
transitions of LLSs are often unsaturated and one may 
(in principal) resolve the H I absorption into multiple 
velocity components, each with a well-measured H I col- 
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uimi density. It would then be possible to explore vari- 
ations in the physical conditions of these components, 
including the gas metallicity. This analysis has not been 
extensively pursued, primarily because of the significant 
observational demands: one requires high S/N, echcllc 
observations covering the full Lyman series. This de- 
mands wavelength coverage A < 4000A for sources at 
z < 3.4 which is observationally challenging. The anal- 
ysis, as well, is both time-consuming and limited by the 
uncertainties of photoionization modeling. With only a 
few exceptions (e.g. Prochaska & Buries 1999), the only 
LLSs analyzed thus far are the rare and special subset of 
LLS that are sufficiently quiescent kinematically to per- 
mit a measurement of the D/H ratio (Buries & Tytler 
1998; Buries ct al. 1999). 

Researchers have focused instead on the integrated 
mctallicities of the so-called super-LLS (SLLS; also 
referred to as sub-DLAs) with TVhi > 10 19 cm~ 2 
(Prochaska 1999; Peroux et al. 2003; Prochaska et al. 
2006). The Lya damping wings for these systems are 
resolved by high-resolution spectroscopy and therefore 
allow a precise constraint on their total Nhi value. In 
turn, the gas metallicity may be estimated subject to 
photoionization corrections. It is rare, however, that the 
full Lyman series is observed or analyzed, and absent the 
TVhi component structure one cannot study metallicity 
gradients. 

This paper marks the second in a series presenting re- 
sults on a survey of z ~ 3 Lyman limit systems. In our 
first paper (O'Meara et al. 2007, ; hereafter Paper I), 
we presented constraints on the H I frequency distribu- 
tion of the SLLSs absorbers with Am = 10 19 cm~ 2 to 
10 20 - 3 cm~ 2 , from a survey of echelle and echellette data. 
In this paper, we present our first detailed analysis of a 
Lyman limit system with Am ~ 10 18 cm~ 2 . Wc demon- 
strate that one can precisely constrain the H I column 
densities of the metal-bearing 'clouds' in such systems 
when provided with high-resolution spectroscopy cover- 
ing the full Lyman series. Furthermore, we measure the 
ionization state, metallicity, and relative abundances of 
the gas through comparisons with models of ionization 
equilibrium. In this fashion, we can also explore metallic- 
ity variation within the LLS, presumably corresponding 
to galactic (sub-Mpc) scales. 

Our target is PKS2000-330, a bright (m v = 17.3) 
quasar at z em — 3.77 (Peterson et al. 1982) discovered by 
Jauncey et al. (1982). Several papers on the subject of 
absorption systems have targeted this quasar previously 
(Hunstead et al. 1986; Murdoch et al. 1986). These au- 
thors characterized the Lya forest and the properties of 
several strong absorption systems along the sightline. In 
their spectra the authors describe four metal-line systems 
with H I absorption indicative of LLSs. Two are a close 
pair at z « 3.19 characterized by typical low- ionization 
lines. Based on a measurement of the O I A1039 line, 
the authors reported that one of these systems displays 
an enrichment level comparable to the Milky Way inter- 
stellar medium (> \ solar). The highest redshift system 
analyzed is the one responsible for the obvious LLS fea- 
ture below ~ 4200A. At z a b s <~ 3.55, this was at the 
time the highest redshift LLS system studied. Limited 
by spectral resolution, the authors presented only a cur- 
sory analysis of its properties. These absorption systems, 



especially the LLS, are the focus of the current paper 4 . 
Throughout the manuscript we assume the solar abun- 
dances reported by Grevesse et al. (2007). 

2. DATA 

We observed PKS2000-330 with the MIKE echelle 
spectrograph (Bernstein et al. 2003) on the Magellan II 
(Clay) Telescope for two 2,400 second exposures on the 
night of September 2, 2004. The MIKE spectrometer 
uses a dichroic centered at A « 5000A to separate the op- 
tical light into two cameras where it is collimated and dis- 
persed onto the SITE 2k x 4k 15/xm detectors. The data 
was reduced using the MIKE Reduction Pipeline 5 (Bern- 
stein et al. 2008) which employs novel techniques for the 
flat fielding, sky subtraction, spectral extraction, wave- 
length and flux calibration of echelle spectroscopy. The 
reduction algorithms also estimate and subtract scat- 
tered light detected between the echelle orders. This 
signal smoothly varies across the image and has a neg- 
ligible contribution compared to the sky, even at blue 
wavelengths. 

The resolution of this data is FWHM~ 10 km s _1 for 
the blue data (3200 - 5000 A) and - 12 km s" 1 for the 
red data (5000 — 9000 A) as measured from a series of 
well-isolated arc lines. The signal-to-noise (S/N) of the 
data is w 25 per 3km s _1 and 4.2km s _1 pixel at 
A =4400A and 6000A, on the blue and red sides respec- 
tively. Because of the differing spectral resolution, wc 
did not coadd the blue and red-side data in the region 
of spectral overlap. Instead, we restrict the analysis of 
features at A < 5000A to the blue-side spectrum and 
A > 5000A to the red-side data. 

The calibrated spectra are presented in Figure 1. It 
is our experience that the relative flux calibration is ac- 
curate to within 15% (Bernstein ct al. 2008). To an- 
alyze the absorption lines in a quasar spectrum, one 
must first estimate its intrinsic continuum. We fit a se- 
ries of b-splines and low-order polynomials to our spec- 
trum of PKS2000— 330 by interactively selecting regions 
where a visible inspection suggests minimal absorption, 
i.e. the regions were assumed to reflect the quasar contin- 
uum. This process is augmented with corrections made 
by hand, generally to maintain a smooth curve. The 
continuum estimate is the leading source of systematic 
error, especially for analysis within the Lya forest. At 
this high redshift, there are few (if any) narrow bands 
of unabsorbed quasar light blueward of the Lya emis- 
sion peak. For spectra with S/N typical of our dataset, 
it is difficult to avoid a systematic underestimate of the 
continuum at these wavelengths. This is especially true 
for data blueward of the rest-frame Lyman limit for the 
quasar (A < 4350A) where absorption from all members 
of the Lyman series may contribute. In Figure 1 we show 
our estimate of the quasar continuum as a light dotted 
line. The reduced and normalized 1-d spectra are avail- 
able at http://www.ucolick.org/^xavier/LLS. 

3. IONIC COLUMN DENSITIES 

In this section, we discuss the techniques and measure- 
ment of ionic column densities for the gas comprising the 

4 The fourth system (at 2 = 3.33) noted by Hunstead ct al. 
(1986) has a poorly constrained TVhi value and we have chosen not 
to analyze its properties in detail. 

5 http://web.mit.edu/~burles/www/MIKE/ 
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Fig. 1.— Magellan/MIKE spectrum of the quasar PKS2000-330 observed on the night of September 2, 2004 (4800s total exposure). 
The upper panel shows the data acquired through the red camera at a dispersion of FWHMsj 12 km S . Strong Lya absorption from 
a complex of super-LLS is indicated at A ss 5100A. The lower panel plots the data acquired through the blue camera at a dispersion of 
FWHM 10 km s — 1 . We mark the Lyman limit system at A fa 4100A corresponding to z ss 3.55. The data have been smoothed in each 
panel by 3 pixels for presentation purposes only (finer detail cannot be resolved by this figure). The dotted line in each panel shows our 
estimate of the quasar continuum which is used to normalize the spectra for absorption-line analysis. 



LLS at z « 3.55 toward PKS2000-330. An analysis of 
the column densities for the SLLSs at z « 3.19 is pre- 
sented in Appendix B. 

3.1. Techniques 

The standard method of characterizing an absorption 
line is with a Voigt profile described by three physi- 
cal parameters: the redshift z, the column density N, 
and the Doppler parameter (b = ay/2). These values 
bear upon the physical quantities that determine the gas 
ionization, most importantly the spatial density, kine- 
matics, and temperature. To fit for these parameters, 
we employed the VPFIT software package kindly pro- 
vided by R. Carswell and J. Webb 6 . For this software 
package, the user supplies (1) spectral regions for analy- 
sis, (2) a list of components with ions specified, (3) ini- 
tial guesses for the parameter values, and (4) the spec- 
tral line-profile parameterized here as a Gaussian with 
FWHM estimated from the ThAr data. In our analysis, 
we have generally avoided spectral regions that are sig- 
nificantly blended with coincident Lya forest absorption 

6 See http://www.ast.cam.ac.uk/~rfc/vpfit.html 



or coincident metal-line absorption. We typically have 
constrained the low-ion 7 components of a system to have 
identical redshift so that the redshift for each set of low- 
ion transitions of a given component are parameterized 
by a single value. This constraint is observationally mo- 
tivated by the good alignment between the absorption 
lines in velocity space. For H I Lyman lines, which suf- 
fer from significant line-blending, we have anchored their 
redshift to the value derived from an independent fit of 
the metal-line transitions. This presumes that the H I ab- 
sorption traces the low-ion profiles in velocity space, but 
the relative column densities are allowed to vary freely 
from component to component. Unfortunately, in the 
LLS presented here the line-blending is too severe to test 
the former assumption empirically. Future analysis of 
the other LLS in our survey will address this point al- 
though we note that such analysis is best performed in 

7 Low-ions are defined to be the first ionization state for a gas 
that has an ionization potential exceeding 1 Ryd (e.g. O , Si+, 
Fe+). These atoms and ions are the dominant ionization state in 
a neutral hydrogen gas where a significant far-UV (hv < 13 eV) 
radiation field is present. 
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the low redshift universe where line-blending with the 
Lya forest, is minimized (e.g. Lehner et al. 2009). For 
lines analyzed within the Lya forest we allow for a sys- 
tematic uncertainty in the continuum placement of 10%. 
This error does not dominate in any of the analysis pre- 
sented in this paper. In general, our fits adopt the fewest 
components which yield reduced y 2 values near unity. 

We complement this line-profile analysis with ionic 
column densities measured from the apparent optical 
depth method (AODM; Savage & Sembach 1996). In the 
AODM, one converts the observed normalized flux pro- 
file (/) to an apparent optical depth array r a = — ln(l//) 
and then sums over the velocity interval of the line-profile 
to calculate the integrated column density. This non- 
parametric technique provides accurate results for well- 
resolved profiles. It also provides conservative upper and 
lower limits for undetected or saturated lines. Through- 
out this paper, we adopt these limits where appropriate. 
Upper limits correspond to 2a statistical uncertainties 
and lower limits correspond to the column densities by 
demanding / > la in every pixel. Finally, the AODM 
analysis is useful for identifying unresolved components 
(hidden saturation) by comparing the column densities 
of a series of transitions with different oscillator strengths 
from a single ion. This analysis is especially appropriate 
for our dataset because the spectral resolution is lower 
than the echelle data commonly used in quasar absorp- 
tion line studies. For the atomic data, we have relied on 
the compilation of atomic data given by Morton (2003). 

3.2. The z » 3.55 LLS Intervening PKS2000-330 

The absorption system at z w 3.55 gives rise to 
the strong Lyman limit observed in the spectrum of 
PKS2000-330 (Figure 1) and, therefore, is the only LLS 
showing the complete Lyman series. As demonstrated 
below, spectral coverage that includes the entire Lyman 
series is crucial to precisely constraining the H I column 
density. We place a lower limit on the total H I col- 
umn density of this LLS by estimating an upper limit to 
the normalized flux shortward of the Lyman limit. In the 
spectral window A = 4000 -4050 A, we measure an upper 
limit to the normalized flux Jll < 0.004 (95% c.l.) where 
we have conservatively allowed for a 50% lower contin- 
uum then presented in Figure 1. This implies a lower 
limit to the total H I column density of the LLS, Nm > 
tll/v{E) = 10 177 cm~ 2 where t ll = ln(l// LL ) and 
a(E) is evaluated at E = 4025A/[911.76A(1 + z LL )] Ryd. 
This assumes that the contribution of coincident absorp- 
tion by the IGM and metal-lines is minor, which is rea- 
sonable for wavelengths longward of the next lower red- 
shift LLSs. We have also ruled out coincident, strong 
Lya absorption at these wavelengths by searching for 
corresponding metal- line absorption (e.g. C IV). 

We can set a conservative upper limit on the total H I 
column density through analysis of the Lya transition. 
Figure 2 presents the spectrum of covering the Lya tran- 
sition for the z ss 3.55 LLS. The dotted line at relative 
flux / ss 3.44 traces the continuum estimated from re- 
gions of minimal absorption within and nearby this spec- 
tral region. Overplotted on the data is the set of Lya 
models corresponding to Nm = 10 181 to 10 18 - 85 cm~ 2 for 
a single absorber centered at z = 3.5523, where the ma- 
jority of H I gas is expected to lie (see below). The profile 
assuming Nm = 10 18 85 cm -2 significantly underpredicts 
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Fig. 2. — Spectral region encompassing the Lya transition of 
the LLS at z pb 3.55. The dotted line at relative flux / m 3.44 
traces the continuum estimated from assumed regions with minimal 
absorption within and nearby this spectral region. Overplotted 
on the data is the set of Lya models corresponding to Nm = 
10 18 1 to 10 18 - 85 cm (solid curve) for a single absorber centered 
at z = 3.5523, where the majority of H I gas is expected to lie in 
this LLS (see below). The profile assuming Nm = 10 18 ' 85 cm - 2 
significantly underpredicts the flux at A rj 5538. 5A. Given that 
the flux at these wavelengths is further depressed by the 'clouds' 
centered at A fa 5537A and 5540A, we adopt a conservative upper 
limit to the Nm value at z = 3.5523 of 10 f8 - 85 cm -2 . 

the flux at A w 5538. 5A, especially if one accounts for the 
additional Lya absorption from the clouds at A w 5537A 
and 5540A. We adopt a conservative upper limit for the 
gas at z = 3.5523 of JV m < 10 18 - 85 cm" 2 . 

A line-profile analysis of the full Lyman series will, in 
principle, give a more precise constraint on the H I col- 
umn density of the LLS. Furthermore, such an analysis 
may also resolve the H I system into multiple velocity 
components. There are several challenges, however, to 
this analysis: (i) the various components of the LLS may 
blend with one another, especially in the higher-order 
Lyman series transitions; (ii) the Lyman series of the 
LLS is severely blended with the lower redshift IGM; 
and (iii) the majority of lines lie on the flat portion of 
the curve-of-growth implying a significant degeneracy be- 
tween b and Nm- We address the first concern (in part) 
by demanding that the strongest H I clouds have iden- 
tical redshift as the strongest low-ion absorption. This 
requirement follows from the assumption that low-ions 
primarily arise in regions of large surface density of H I 
gas that can self-shield the material from ionizing pho- 
tons. We address the second point by avoiding severely 
blended lines. On the third point, we consider physically 
plausible values for the Doppler parameter and adopt the 
corresponding TVhi constraints. 

In Figure 3, we present a series of metal-line transi- 
tions for the LLS at z « 3.55. We have divided the line- 
profiles into three subsystems: (A) —230 km s _1 < Sv < 
-70 km s" 1 , (B) -70 km s" 1 < Sv < 70 km s -1 , and (C) 
70 km s" 1 < Sv < 210 km s -1 relative to z LLS = 3.550. 
The velocity interval of the subsystems was chosen to 
separate the two significant low-ion complexes from each 
other (B and C) and also to separate strong C IV ab- 
sorption (A) from the other subsystems. Subsystems B 
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Fig. 3. — Velocity plot of the principal metal-line transitions detected for the LLS at z ft: 3.55 (v = km s — 1 corresponds to z = 3.550). 
In each panel, the vertical (green) dashed lines designate three subsystems (A,B,C) that we have designated. Blends with coincident 
absorption (e.g. Lya lines) are denoted as (orange) dashed spectra. Overplotted on the data is the best-fit solution from a line-profile 
analysis of the data using the VPFIT software package. Arrows in the top panels locate the centroids of the velocity components (Table 1). 



and C show significant low-ion absorption and we ex- 
pect that these contain the majority of the H I gas. We 
have independently fit the metal-line profiles of the var- 
ious subsystems using the VPFIT software package. We 
have forced the components comprising each subsystem 
to have identical redshift for different low-ions but have 
allowed the Doppler parameter values to vary because the 
line widths of the high-ion transitions appear systemati- 
cally wider than those of the low-ions. The best-fit solu- 
tions and la error estimates are presented in Table 1 and 
the models are overplotted on the data in Figure 3. The 
low and high-ion absorption in each subsystem is well 



modeled by components having identical velocity, but the 
analysis also suggests that the high-ion gas has system- 
atically larger Doppler parameters (b > 20 km s _1 ) than 
the low-ion gas (6 < 15 km s _1 ). This conclusion, how- 
ever, is sensitive to the details of the component struc- 
ture assumed in the analysis. It is possible to achieve 
an acceptable model where the 6-values are identical be- 
tween the low and high-ions, but this would require the 
inclusions of yet further components. The profile model 
presented here is the simplest one that reproduces the 
observations. We return to this point when discussing 
photoionization models of the system (§4). Finally, we 
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TABLE 1 

VPFIT SOLUTIONS FOR THE LLS AT z = 3.550 



Subsystem 


z 


a(z) 

(io- 5 ) 


v a 

(km s -1 ) 


Ion 


b 

(km s -1 ) 


a(b) 
(km s 1 ) 


log N 


a(N) 


A 


3.547932 


0.5 


-136 


H I 


17.45 


3.27 


14.88 


0.10 










C IV 


14.09 


0.60 


13.95 


0.03 










Si IV 


10.90 


1.34 


13.24 


0.04 










Si III 


16.68 


1.32 


13.00 


0.03 


B 


3.549521 


1.1 


-32 


C IV 


3.00 


1.56 


13.23 


0.16 










Si IV 


3.00 


9.69 


12.09 


0.17 


B 


3.549938 


0.5 


-4 


H I 


19.38 


4.53 


16.82 


0.20 










C IV 


27.48 


1.53 


13.89 


0.03 










CII 


10.01 


0.64 


14.09 


0.04 










OI 


9.24 


3.78 


13.34 


0.09 










Al III 


14.10 


8.28 


12.11 


0.17 










Si II 


8.97 


1.50 


13.27 


0.04 










Si IV 


17.97 


1.53 


13.11 


0.03 


C 


3.550751 


3.8 


+50 


CIV 


17.49 


4.60 


13.13 


0.09 


C 


3.552304 


0.4 


+152 


HI 


25.00 


2.50 


18.00 


0.25 










CIV 


23.71 


1.11 


13.65 


0.03 










CII 


12.86 


0.42 


14.00 


0.03 










Al III 


8.58 


0.42 


11.96 


0.18 










Si IV 


18.67 


0.48 


13.64 


0.03 










Si II 


8.41 


0.42 


12.84 


0.11 



a Velocity relative to 2 = 3.550. 

note that the regions containing the transitions of the 
O VI and N V doublets arc hopelessly blended with co- 
incident IGM features. 

The redshifts of the strongest low-ion component in 
each of the subsystems were then taken as an input con- 
straint for the H I Lyman series analysis. In Figure 4, 
we present a subset of the Lyman series transitions and 
also a portion of the data near the Lyman limit. It is 
evident that the strongest, highly saturated members of 
the Lyman series offer only a weak constraint on the H I 
column densities of the subsystems. We first visually 
identified spectral regions that were most likely to con- 
strain the strongest H I components of subsystems A, B, 
and C. We then inputted to VPFIT these regions, a set 
of H I absorption lines at the redshifts of each subsystem 
[za = 3.547932, z B = 3.549938, z c = 3.552304), and ad- 
ditional H I lines within subsystem B (with zc > z > Zb) 
to model additional absorption evident in the Lyman se- 
ries. The non-zero flux at Sv w 50 km s -1 in the high- 
order Lyman scries lines (e.g. H I 918, 919) indicates 
that the integrated column densities of these additional 
components is significantly lower than the H I column 
density of subsystem C and possibly subsystem B. We 
can set a conservative upper limit to the total Am of 
this gas by integrating the apparent optical depth of the 
H I 918 profile in the interval Sv = +17 to +80 km s~\ 
finding A ffl < 10 16 - 9 cuT 2 (95% c.L). This value is to 
be considered an upper limit because there may be sig- 
nificant line-blending with lines from lower redshift IGM 
systems. 

Our trials with the VPFIT package yielded a set of 
models with acceptable reduced x 2 values, but we found 
that the parameters of the additional components in sub- 
system B are very poorly constrained and that we could 
not achieve a unique solution for the full set of profiles. 
Therefore, we used the range of solutions to guide a 'by- 
eye' fitting analysis of the components at z — zb and zc- 
The parameters for the gas at z = zb are well constrained 
by the higher order Lyman series (H I 918,919,920,923): 



logiV H i = 16.8 ±0.2, b m = 19.4 ±4.5 km s" 1 . There 
is a weak dependence of these values on our assumed 
parameters for absorption at Sv w +20kms _1 which 
contributes to our error estimate. Nevertheless, the to- 
tal JVhi value for the gas in subsystem B is much less 
than that required to explain the observed absorption at 
A < 4150A. We conclude that the majority of H I gas in 
the LLS at z w 3.55 is associated with subsystem C. 

The detection of positive flux at Sv w +180 km s -1 in 
H I 937 and at Sv w +80 km s -1 in H I 918,919 sets an 
upper limit to the combined Ahi and b m values of an 
H I component centered at z — zc- Formally, the data 
permit Am values as large as 10 18 85 cm _2 provided b- 
values less than 10km s _1 . We adopt a prior of b m > 
20 km s _1 , however, based on the following physical ar- 
guments. First, we demonstrate in the following section 
that the gas related to subsystem C is highly ionized. 
Even if we assume a soft (i.e. stellar) ionizing radiation 
field, the implied gas temperature is T > 15, 000K giving 
b > 16 km s _1 . Second, systems with b < 20 km s -1 are 
extremely rare, as evidenced by the rarity of D/H mea- 
surements from the IGM. Third, the C II transitions have 
Dopplcr parameters b > lOkms -1 which set a lower 
limit to the Doppler parameter of the H I gas. Finally, 
a model with b < 20 km s -1 would underpredict the ab- 
sorption at Sv w — 45kms~ 1 and +60 km s -1 for the 
entire Lyman series. Therefore, one would need to intro- 
duce an additional 10 or more absorption lines (including 
many at unrelated redshifts) to reproduce the observa- 
tions. Together, these points motivate a lower limit to 
6gj of w 20 kms- 1 . 

With this prior on b m , the data require Am < 
18.5 cm -2 and we adopt a best-estimate of Nm — 
10 i8.o±o.25 cm -2 and 6 c = 25 ± 2.5 km s -1 for subsys- 
tem C. We caution that subsequent analysis should not 
treat the Am or ^-values for subsystem C as following 
a normal distribution with the reported uncertainties. 
Instead, we recommend that the central value be consid- 
ered uniformly distributed. But we also note that the 
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Fig. 4. — Velocity plots of the Lyman limit (upper-left panel) and a set of H I Lyman series transitions for the LLS at z ~ 3.55 toward 
PKS2000— 330. Arrows designate the centroids of the various velocity components used in the analysis (Table 1). The solid blue curve 
shows the convolved line-profile for each Lyman line individually, while the dotted gray curve plots the convolved solution for all Lyman 
lines (inspect the panel showing the H I 918 transition where multiple Lyman lines overlap). The vertical dashed lines in the panels showing 
H I 915b and 930 designate subsystems A, B, and C. The velocity v = km s — 1 corresponds to an arbitrary z = 3.550. 



TVhi value cannot be significantly less than 10 cm 
because of the combined constraints from the Lyman 
limit and subsystem B. 

The full set of constraints on the Nm values for the 
subsystems comprising this LLS is summarized in Ta- 
ble 2. Our favored solution is overplotted on the data 
in Figure 4. It is evident that the strongest members of 
the Lyman series (a, (3, 7, 6) offer little constraint on the 
A^hi values of subsystems B and C. The best constraints 
come from higher order lines; these are presented in Fig- 
ure 5. In this figure, we present the favored solution 
(middle panel) and 2a departures in the iVui and b val- 



TABLE 2 

CONSTRAINTS ON JV H i FOR THE LLS AT z 



3.550 



Subsystem 


Feature 


log AT HI 


Assumptions 


A+B+C 


Lyman Limit 


> 17.7 


Minimal IGM blending 


C 


Lya 


< 18.85 


None 


A 


Ly/5, Ly7, Ly<5 


15.0 ±0.1 


Minimal IGM blending 


B 


Ly8-ll 


16.8 ±0.2 


Minimal IGM blending 


C 


Ly5, Lyl0-ll 


18.0 ±0.25 


b > 20 km s" 1 



ues (incremented in opposition to minimize the change to 
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TABLE 3 

IONIC COLUMN DENSITIES FOR THE LLS AT z = 3.55 

Ion A rcst log/ Vf nt logJVAODM log iVyPFIT log iV a dopt 

(A) (km s" 1 ) 



SUBSYSTEM A 



HI 


1215 


6701 


-0.3805 








14.88 ± 0.20 


CII 


1334 


5323 


-0.8935 


-186, -86 


< 12.75 




< 12.75 


CIV 


1548 


1950 


-0.7194 


-186, -96 


13.91 ±0.03 


13.95 ±0.03 


13.95 ± 0.03 




1550 


7700 


-1.0213 


-186, -96 


13.93 ±0.03 






OI 


1302 


1685 


-1.3110 


-166, -96 


< 13.00 




< 13.00 


Al III 


1854 


7164 


-0.2684 


-186, -86 


< 12.12 




< 12.12 




1862 


7895 


-0.5719 


-186, -86 


< 12.37 






Si II 


1304 


3702 


-1.0269 


-186, -86 


< 12.80 




< 12.72 




1526 


7066 


-0.8962 


-166, -96 


< 12.72 






Si III 










13.00 ±0.03 


13.00 ± 0.03 


Si IV 


1393 


7550 


-0.2774 


-186, -96 


13.14 ±0.03 


13.24 ±0.04 


13.24 ±0.04 




1402 


7700 


-0.5817 


-186, -96 


13.23 ±0.03 






Fell 


1608 


4511 


-1.2366 


-186, -86 


< 12.93 




< 12.93 


SUBSYSTEM B 


HI 


1215 


6701 


-0.3805 








16.82 ± 0.20 


CII 


1036 


3367 


-0.9097 


[-47, 33] 


14.12 ±0.03 


14.09 ±0.04 


14.09 ± 0.04 




1334 


5323 


-0.8935 


[-57, 43] 


14.03 ± 0.03 








1335 


7077 


-0.9397 


[-47, 23] 


< 13.36 






CIV 


1548 


1950 


-0.7194 


[-57, 33] 


13.89 ±0.03 


14.03 ±0.03 


14.03 ± 0.03 




1550 


7700 


-1.0213 


[-57, 33] 


14.05 ± 0.03 






OI 


1302 


1685 


-1.3110 


[-37, 33] 


13.35 ± 0.07 


13.34 ±0.09 


13.34 ±0.09 


Al III 


1854 


7164 


-0.2684 


[-57, 43] 


12.19 ±0.13 


12.11 ±0.17 


12.11 ± 0.17 




1862 


7895 


-0.5719 


[-57, 43] 


< 12.38 






Si II 


1304 


3702 


-1.0269 


[-57, 43] 


13.21 ±0.06 


13.27 ±0.04 


13.27 ±0.04 




1526 


7066 


-0.8962 


[-37, 33] 


13.29 ±0.03 






Si IV 


1393 


7550 


-0.2774 


[-57, 33] 


13.14 ±0.03 


13.15 ±0.03 


13.15 ± 0.03 




1402 


7700 


-0.5817 


[-57, 33] 


13.18 ±0.03 






Fell 


1608 


4511 


-1.2366 


[-57, 43] 


< 13.09 




< 13.09 


SUBSYSTEM C 


HI 


1215 


6701 


-0.3805 








18.00 ±0.25 


CII 


1036 


3367 


-0.9097 


[112, 192] 


14.20 ±0.03 


14.00 ±0.03 


14.00 ±0.03 




1334 


5323 


-0.8935 


[102, 202] 


13.97 ±0.03 








1335 


7077 


-0.9397 


[122, 192] 


12.90 ±0.10 






CIV 


1548 


1950 


-0.7194 


[102, 192] 


13.65 ± 0.03 


13.65 ±0.03 


13.65 ± 0.03 




1550 


7700 


-1.0213 


[102, 192] 


13.64 ±0.03 






OI 


1302 


1685 


-1.3110 


[122, 192] 


< 13.01 




< 13.01 


Al III 


1854 


7164 


-0.2684 


[102, 202] 


< 12.13 


11.96 ±0.18 


11.96 ± 0.18 




1862 


7895 


-0.5719 


[102, 202] 


< 14.04 






Si II 


1304 


3702 


-1.0269 


[112, 192] 


12.79 ±0.14 


12.84 ±0.11 


12.84 ±0.11 


Si III 


1206 


5000 


0.2201 


[102, 202] 


> 13.59 




> 13.59 


Si IV 


1393 


7550 


-0.2774 


[102, 192] 


13.60 ±0.03 


13.64 ±0.03 


13.64 ±0.03 




1402 


7700 


-0.5817 


[102, 192] 


13.62 ±0.03 






Fe II 


1608 


4511 


-1.2366 


[102, 202] 


< 12.91 




< 12.91 



a Velocity interval for the AODM relative to z = 3.550000. 



the model). The shaded regions identify the pixels that 
most constrain these parameters. A few points should 
be emphasized. First, the left-hand panels correspond 
to a model with A%i = 10 17 ' 5 cm~ 2 which is ruled out 
by observations of the Lyman limit (see above). Second, 
the model in the right-hand panel significantly underpre- 
dicts the absorption at Sv s=a +65 km s _1 . Full coverage 
of the Lyman series and the Lyman continuum region 
at a relatively high S/N ratio has constrained the Nm 
value of this LLS. Such analysis demands echelle spectra 
of the full Lyman series can minimize blending with the 
z <~ 3 IGM and the foreground LLSs. The full set of ionic 
column densities for the three subsystems is provided in 
Table 3. 

Before concluding our discussion on the H I absorp- 
tion of the LLS at z i=a 3.55, we remark that our model 



is not a good fit to the spectrum at A w 4160 — 4175A. 
In this spectral range, which corresponds to the Ly-12, 
Ly-13, and Ly-14 transitions of the LLS, our model un- 
derpredicts the absorption at velocities in between the 
main components (i.e. away from the line centers). This 
implies either an additional source of nearly continuous 
opacity unrelated to the LLS at z — 3.55 or that we have 
significantly overestimated the quasar continuum flux at 
these wavelengths. Comparing to the published spec- 
trum of Barthcl ct al. (1990), we note a similar drop in 
the quasar flux at these wavelengths which argues against 
the latter explanation. Our expectation, therefore, is 
that the majority of additional opacity is from unidenti- 
fied H I lines from the intergalactic medium at z < 3.55. 
In any case, the measurements of TVhi for this LLS are 
well constrained by the Lyman lines redward of 4180A. 
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Fig. 5. — Zoom-in on the Lyman series lines that most tightly constrain the column densities and b- values for the H I components 
comprising subsystem C. The center panels show the favored model (solid blue curve) line for the central values of our line-profile solution 
while the left/right panels show 2<r departures from the central value chosen to minimize changes in the model, i.e. an increase in JVhi is 
offset by a decrease in the b-value. The shaded boxes show the portion of the data that has greatest constraint on the analysis. 



4. PHOTOIONIZATION MODELING 

At z > 2, gas clouds with H I column densities of less 
than 10 19 cm -2 are predicted to be photoionized by the 
extragalactic ultra-violet background (EUVB) radiation 
field unless one assumes an extraordinarily high volume 
density (e.g. Viegas 1995). In this photoionized gas, the 
ions observed (especially the low-ions and H I) may rep- 
resent only a trace quantity of the total gas present. To 
determine the gas abundance and also assess the physical 
conditions, it is necessary to model the ionization state of 
the gas. We have performed this modeling for the LLSs 
along the sightline to PKS2000-330 using the Cloudy 
software package (Ferland et al. 1998). In practice, we 



have calculated ionization models for a series of plane- 
parallel slabs with constant gas density ur and with JVhi 
constrained to match the central values derived from our 
line-profile analyses. We have adopted the EUVB radia- 
tion field calculated by Haardt & Madau (in prep.; a.k.a. 
CUBA) assuming contributions from both quasars and 
UV-bright galaxies. We run models with a wide range 
of ionization parameters U = $/rz.jjc, which is akin to 
adopting an intensity $ for the EUVB 8 and varying the 
Hydrogen volume density nn- Finally, we compare ob- 

8 This intensity has been poorly constrained by empirical obser- 
vation, but see Dall'Aglio et al. (2008). 
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served ionic ratios against the model predictions to con- 
strain the ionization parameter U under the assumptions 
of our simplistic photoionization model. Ideally, we con- 
sider only pairs of ions from the same element to min- 
imize dependence on the assumed intrinsic abundances. 
But we also find valuable constraints from ions of dif- 
ferent elements, even allowing for significant departures 
from non-solar relative abundances. We caution that the 
single-phase models considered here are overly simplistic. 
In Appendix A, we present two-phase models to explore 
the impact of more complex scenarios. 

As noted above, an absorption system with iVni -C 
10 20 cm -2 is predicted to be highly ionized at z > 3 by 
the EUVB radiation field alone. One signature of signif- 
icant photoionization is the presence of strong high-ion 
absorption (C IV, Si IV) relative to low-ion species. In- 
deed, we find that the ratios of C+/C+ 3 and Si+/Si+ 3 
are of order unity for subsystems B and C and much less 
than one for subsystem A. Under the assumption of a 
single-phase model, these observations indicate the gas 
in this LLS is predominantly ionized. As noted in § 3.2, 
however, the high-ion gas may have systematically larger 
Dopplcr parameters than the low-ion species. This might 
indicate that the gas occurs in a distinct region, physi- 
cally separated by a large distance from the low-ion gas. 
We consider this hypothesis improbable because (i) it 
would require a remarkable coincidence for the low-ions 
and high-ions to arise in truly distinct 'clouds' yet share 
the same velocity; and (ii) the low TVhi of this LLS re- 
quires high ionization fractions and corresponding Si IV, 
C IV absorption unless one invokes a very large gas den- 
sity tih > 0.1 cm~ 3 (see Appendix A). We attribute 
the differences in 6-values to the effects of photoioniza- 
tion from external sources onto an optically thick, self- 
shielding gas. Optical depth effects imply that the lower 
ionization states lie preferentially toward the middle of 
the absorbing gas while the high-ion gas arises in the 
outer regions. In this scenario, one expects a tempera- 
ture gradient which would imply higher 6-values for the 
high-ions. One can also allow for (and may even ex- 
pect) a gradient in velocity dispersion that would further 
boost the Doppler parameters. A full treatment of such a 
model is beyond the scope of this paper. We proceed as- 
suming a single-phase (constant density) , optically thick 
model and restrict the ionic ratios to the gas at the veloc- 
ities Sv = -162,-30, +126 km s" 1 as listed in Table 1. 
In Appendix A, we explore how the results vary if we 
assume that the majority of the high-ion absorption is 
unrelated to the low-ion gas. We find that the principal 
scientific conclusions are qualitatively the same. 

In Figure 6 we present the predicted ratios for pairs 
of ions observed in the LLS. Overplotted on the curves 
are the observed values (assuming a minimum uncer- 
tainty of 0.15 dex) where the horizontal error bars in- 
dicate the implied constraints for the ionization param- 
eter. Allowing for the significant uncertainties inherent 
to this modeling (simplistic geometry, uncertain atomic 
data, etc.), the observations are broadly consistent by 
taking \ogU A = -1.8 ± 0.3, log U B = -2.8 ± 0.2, and 
logt/c = —2.4 ±0.2 for the three subsystems. These val- 
ues are typical of the ionization parameters derived for 
a handful other LLSs to date (e.g. Buries & Tytler 1998; 
Prochaska 1999; Prochaska & Buries 1999; Prochaska & 
Hennawi 2009). These estimated ler uncertainties are 
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Fig. 6. — These curves show the predictions from a series of 
Cloudy photoionization models for the column density ratios of Si, 
C and O atoms and ions. The models are parameterized by the 
ionization parameter U which corresponds to the ratio of ionizing 
photons to hydrogen atoms per unit volume. In each panel, we 
have tuned the Cloudy models to match the observed A^jji values 
and metallicities of subsystems A, B, and C from the LLS at z rj 
3.55. Overplotted on the curves are the observational constraints. 
We have used these to estimate the ionization parameter for each 
subsystem (designated by the shaded regions in each panel). 

dominated by systematic (not statistical) uncertainty. 
Note that the results do not change even if we assume 
that 90% of the high-ions arise in an alternate phase (Ap- 
pendix A). This is because the curves involving Si IV and 
C IV are steep functions of log U. We also note that the 
Al + /Si + observations indicate a significantly ionized gas 
even for our assumed solar intrinsic Al/Si abundance. 

All of the observed ionic ratios are roughly consistent 
with the assumed U values save one: the O /Si + ratio 
for subsystem B. In this subsystem, the observed col- 
umn density ratio is log[A^(O )/A(Si+)] = +0.03 ± 0.1 
yet the photoionization model predicts iV(O ) <§C A(Si + ) 
for solar relative abundances and \ogUs > — 4. At our 
preferred value of log Ub — —2.80, the disparity between 
model and observation exceeds one order of magnitude. 
We have carefully considered the possibility that the ab- 
sorption line at A — 5925A identified as O I 1302 is a 
mis-identification, especially in light of its non-detection 
in subsystem C. The line, however, is not associated with 
the other subsystems, the SLLSs at z ~ 3.19, the Mg II 
absorber at z « 2.03, C IV 1548, nor any other fre- 
quently observed doublet of the IGM. We are confident 
in the identification, therefore, but are challenged to ex- 
plain the offset between observation and model in the 
O /Si + ratio. One possibility is that the gas has a highly 
super-solar O/Si ratio, i.e., [O/Si] ^> 0. This runs con- 
trary, however, to theoretical expectation (and empirical 
observations) that a-elements roughly trace one another 
because they are both produced mainly by massive stars. 
At most, one might allow for departures from the solar 
abundances of ~ 0.3 dex. 
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TABLE 4 
ELEMENTAL ABUNDANCES 
FOR THE LLS AT 2 « 3.55 



Ion [X/H] [X/Si+] 





Subsystem 


A a 


c+ 


< +0.08 


< -0.91 


C +3 


-0.33 ± 0.50 


-1.32 ± 0.63 


o° 


< +4.00 


< +3.01 


A1++ 


< +0.53 


< -0.47 


Si+ 


< +0.99 




Si+ 3 


-0.16 ±0.24 


-1.15 ± 0.37 


Fe+ 


< +3.86 


< +2.87 


Subsystem B 6 


C+ 


-0.66 ±0.20 


-0.05 ±0.20 


C +3 


-0.10 ±0.48 


+0.50 + 0.49 


o° 


> -0.14 


> -1.08 


A1++ 


-0.95 ±0.20 


-0.35 ± 0.20 


Si+ 


-0.61 ±0.20 




Si+ 3 


-0.73 ± 0.40 


-0.13 ± 0.40 


Fe+ 


< +0.45 


< +1.05 


Subsystem C c 


C+ 


-1.89 ±0.20 


+0.30 + 0.20 


C +3 


-2.34 ±0.43 


-0.14 ±0.35 


o° 


< -0.36 


< +1.83 


A1++ 


-2.26 ±0.20 


-0.07 + 0.20 


Si+ 


-2.19 ±0.20 




Si+ 3 


-1.85 ±0.33 


+0.34 + 0.25 


Fe+ 


< -0.65 


< +1.54 



Note. — In all cases, we have as- 
sumed a minimum error of 0.2 dex due 
to systematic errors in the photoion- 
ization modeling. 

a Assumes a Cloudy photoionization 
model with N m = 10 16 ' 7 cm -2 , 
[M/H] = -0.50 and log U = -1.8 ± 
0.3 dex 

b Assumes a Cloudy photoionization 
model with N m = 10 16 ' 9 cm -2 , 
[M/H] = -0.50 and log U = -2.8 ± 
0.2 dex 

c Assumes a Cloudy photoionization 
model with Ahi = 10 17 ' 9 cm -2 , 
[M/H] = -1.50 and log U = -2.4 ± 
0.2 dex 

The predicted low value for the 0°/Si + ratio from 
our models with \ogU > —4 stem primarily from the 
large cross-section of 0° to extreme ultraviolet and x- 
ray photons (Sofia & Jenkins 1998). If one adopts a 
softer spectrum (i.e. one absent the influence of quasars), 
the 0°/Si + ratio tends toward the intrinsic O/Si abun- 
dance, e.g. log[V(0°)/V(Si+)] w +1 for solar rela- 
tive abundances. It is possible that the observed ratio, 
log[V(0°)/V(Si+)] w 0, indicates a softer radiation field 
than we have adopted. This would require an intense 
and local source of radiation, e.g. the UV flux from a 
star-burst galaxy. Deep imaging of the field does not re- 
veal any nearby bright source of UV flux (O'Meara et al. 
2006). We proceed by adopting the 'quasar+galaxy' 
EUVB model and report the O/H abundance as a lower 
limit from the observed 0°/H° ratio. Interestingly, this 
yields [O/H] > —0.1 dex, a limit that lies three times 
above the ionization corrected [Si/H] value. We note that 
a softer ionizing spectrum would also imply a higher Si/H 



TABLE 5 

SUMMARY OF PROPERTIES FOR THE LLS AT 2 3.55 



Property 


A 


B 


C 


log (AW cm" 2 ) 


14.9 + 0.20 


16.82 ±0.20 


18.0 ±0.25 


log U 


-1.8 ± 0.3 


-2.8 ±0.2 


-2.4 ±0.2 


log(l - x) a 


_o 4 o+0.36 
•5-^-0.35 


-2 46+ ' 21 


c\ f?Q~l~0.23 
— z. (O_ 22 


log(A H /cm- 2 ) 


18.30 ± 0.41 


19.28 ±0.29 


20.73 ± 0.34 


n e (cm -3 ) 




< 13.0 


< 5.0 


n H (cm -3 ) 




< 13.0 


< 5.0 


<(pc) 




> 0.5 


> 34.8 


[O/H] 


< 1.46 


> -0.2 


< -1.65 


[Si/H] 


-0.10 ±0.35 


-0.61 ±0.20 


-2.19 ± 0.29 


[C/H] 


-0.17 ±0.60 


-0.66 ±0.20 


-1.89 ±0.25 


[Fe/H] 


< 3.78 


< 0.45 


< -0.65 



Note. — Chemical abundances [X/H] assume the photoion- 
ization models as described in Table 4 except for O/H which 
adopts no correction. 

a The ionization fraction x is defined as H+/H. 



abundance for subsystem B. Table 4 presents the abso- 
lute and relative abundances for the LLS adopting the 
ionization corrections derived from the 'quasar-galaxy' 
EUVB photoionization model with the exception of oxy- 
gen in subsystem B where we adopt conservative lower 
limits based on the observed 0°/H° ratio. We discuss 
these results at further length in the following section. 

Under the assumption of detailed balance, we can set 
an upper limit to the electron density of the gas from the 
non-detection of C II* 1335 absorption in subsystems B 
and C. Adopting an electron temperature T e = 25, 000K, 
which is appropriate for this photoionized gas, we have 
n e = 117/[2V(C+ =1/2 )/V(C+ = 3 /2 )-l]. The upper limits 

to V(C + =3/2 ) from the non-detections of the C II* 1335 

transition imply n e < 13 cm -3 and < 5 cm -3 for sub- 
systems B and C respectively 9 . Because the gas is pre- 
dominantly ionized, we infer the same upper limits to 
the hydrogen volume densities njj. These results are 
independent of the assumed ionization model. Finally, 
we can estimate lower limits to the characteristic sizes 
I = Nh/tih of the 'clouds' comprising subsystems B and 
C: £b > 0.5 pc and lc > 28 pc, respectively. 

5. DISCUSSION 

In the previous sections, we presented measurements of 
the ionic column densities of the Lyman limit system at 
z w 3.55 along the sightlinc to PKS2000-330. We then 
compared these values against photoionization models 
of plane-parallel gas slabs to infer the ionization state 
of the gas and thereby estimate physical conditions and 
chemical abundances for the subsystems. A summary 
of the key properties is presented in Table 5, under the 
assumption of a simplistic, single-phase photoionization 
model. Appendix A discusses how these results change 
if we adopt more complex models. 

A principal result of this work is the precise mea- 
surement of the H I column density of a LLS with 
Vhi ~ 10 18 cm -2 through a combined analysis of the 
Lyman limit and the full set of Lyman series lines. We 
have demonstrated that echelle spectra of the full Lyman 
series and the Lyman limit can constrain Vhi in the dom- 

9 Note that LLSs with larger C+ column densities should provide 
tighter (more meaningful) constraints on n e . 
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inant low-ion components (assuming identical redshifts) 
to several tenths dex for an absorber which lies on the flat 
portion of the curve-of-growth and which is significantly 
blended with the z ~ 3 IGM and several foreground 
LLSs. Presently, there is very weak empirical constraint 
on the frequency distribution of iV H i values f(N m ) for 
the LLS and Ahi < 10 19 cm~ 2 . The results presented 
here give confidence that one can establish f(Nm) for 
the LLS provided a large survey of systems with echelle 
observations using the full Lyman series. This forms the 
basis of an observational survey we have recently com- 
pleted. 

More importantly, we have resolved the LLS system 
into several metal-bearing velocity intervals (termed sub- 
systems A, B, and C) and independently constrained 
their A%n values, ionization state, and physical proper- 
ties. This enables a probe of metallicity variations at 
Sv < 100 km s -1 velocity separation and in gas asso- 
ciated with a single galaxy halo or galaxy-scale struc- 
ture 10 . As Table 5 reveals, these subsystems exhibit 
large differences in enrichment (> ldex). Although we 
may dismiss subsystem A from the discussion because 
of its much lower A%t value, the metallicities of subsys- 
tems B and C alone diverge by over an order of mag- 
nitude. Although large dispersions in gas metallicities 
have been reported for the intergalactic medium (e.g. 
Schaye et al. 2003) and individual gas-rich galaxies at 
high z (Prochaska et al. 2003; Peroux et al. 2007), these 
studies refer to gas in structures separated by very large 
distances IMpc). Our results, which mark the first 
detailed analysis of multiple components in a LLS drawn 
from a homogeneous sample 11 , indicate a similar dis- 
persion holds down to galactic scales. This conclusion 
resembles that of Schaye et al. (2007), who suggest the 
IGM is inhomogcneously enriched by a population of high 
metallicity, H I-deficient absorbers 12 . The large metallic- 
ity variations are also consistent with the small physical 
sizes implied for metal-line systems from sightline studies 
of lensed quasars (Rauch et al. 2001). 

We have also analyzed the metallicity of the SLLSs at 
z w 3.19 toward PKS2000-330 (Appendix B, Table B3). 
These are separated by several hundred km s _1 and 
are more likely to correspond to multiple galactic ha- 
los and/or large-scale structures. The observational re- 
sults for the LLS, its subsystems, and the SLLSs are pre- 
sented in Figure 7 where the metallicity values are plot- 
ted against the measured H I column densities. At face 
value, the figure reveals a trend of declining metallicity 
with increasing A?hi value. Although this is an accurate 
account of the systems analyzed in our paper, the trend 
is biased by our having focused on metal-bearing gas. 
At the lower Nm values, especially, there are H I compo- 
nents in the LLS (and presumably the SLLSs) that ex- 
hibit only weak or non-discernible metal-line absorption. 
An obvious example of this is the additional components 

10 The association of a Lyman limit system with a single galactic 
structure is well-motivated by the rarity of LLSs along individual 
sightlines. 

11 There are also reports of metallicity variations in the few 
special LLS that permit D/H analysis (Tytler et al. 1996; Kirkman 
et al. 2003). Prochaska & Buries (1999) also presented similar 
results for two subsystems in the partial LLS at z = 1.92 toward 
J2233— 606 noting a difference of at least 0.5 dex in [C/H]. 

12 One may note that subsystem A could be a higher density 
analog of these H I deficient metal-line systems. 
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Fig. 7. — Silicon metallicity [Si/H] versus H I column density 
Nm as measured for the three subsystems (A,B,C) of the LLS at 
2 ~ 3.55 and the three SLLSs (a,b,c) at z ~ 3.19. For each sys- 
tem, the figure also presents the Afjj-weighted, mean metallicity 
< [Si/H] > as an additional error interval marked with an 'x'. 
There is an apparent decrease in [Si/H] with increasing Ahi value. 
This is partly due to observational bias: metal-poor systems with 
low JVht value were not considered in this paper. The large dis- 
persion in metallicity between LLSs and also in the velocity com- 
ponents comprising a given LLS is a robust result of our analysis. 

within subsystem B of the LLS (Table 1). These have 
estimated Nm values of nearly 10 17 cm -2 yet show es- 
sentially no low-ion absorption and only weak high-ions 
(Figure 3). We have not estimated metallicities for these 
components, in part because we lack constraints on the 
ionization state. Nevertheless, this gas likely has signif- 
icantly lower metallicity than subsystem B and possibly 
even subsystem C. It is our expectation, therefore, that 
the low Njji, low metallicity region of Figure 7 is popu- 
lated by gas related to the LLS. The most robust result 
regarding the metallicities, therefore, is that the gas in 
LLSs exhibits a large dispersion (i.e. greater than 1 dex). 
This holds not only for the subsystems comprising the 
LLS but also among the various LLSs studied here. 

The most natural interpretation of the large metallicity 
variations is the incomplete mixing of metals in the asso- 
ciated galactic-scale structure. Presumably the gas has 
previously cycled through the neutral ISM of a nearby 
galaxy 13 (i.e. a damped Lya system) and subsequently 
was transported to the surrounding region. One plausi- 
ble explanation is that this transport process does not 
efficiently mix the gas. Another possibility is that the 
large metallicity variations represent a homogeneously 
enriched material that has mixed inhomogeneously with 
'fresh' primordial gas. In both scenarios, we conclude 
that mixing is incomplete on galactic scales and that 
large metallicity variations may occur in the galaxy/IGM 
interface. 

The large metallicity variation observed within the 
LLS may also have important implications for larger 
Ahi systems where the H I absorption is not resolved 
and, therefore, metallicity gradients cannot be assessed. 

13 An inspection of deep optical images of the field surround- 
ing PKS2000— 330 does not reveal any obvious z ~ 3.5 galactic 
counterpart (O'Meara et al. 2006). 
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Metallicity variations within the SLLSs and DLAs could 
occur if the sightlines sample a metallicity gradient 
within the galaxy and/or multiple gas phases with vary- 
ing enrichment levels surrounding high z galaxies (Wolfe 
& Prochaska 2000; Fox et al. 2007; Lehner et al. 2008). A 
large metallicity dispersion within DLAs runs contrary, 
however, to the relatively small dispersion observed in 
the relative abundances of a/Fe, N/O, etc. (Prochaska 
2003; Dessauges-Zavadsky et al. 2006). Nevertheless, we 
are motivated to explore metallicity variations in a larger 
sample of LLSs and to seek trends with the total Ahi 
value of the system. This is of special interest given that 
the A H value of the LLS and SLLS is representative of 
DLAs. The principal physical difference could simply be 
the ionization parameter of the gas. 

Figure 7 also presents the An-weighted, mean metal- 
licity integrated over the entire LLS and SLLS systems, 
i.e. 

< [Si/H] >= (1) 

i 

where A H was derived from photoionization modeling. 
Regarding the LLS at z ~ 3.55, we find that its mean 
metallicity is dominated by subsystem C which has the 
lowest level of enrichment but the highest Ah value 14 . 
If we had derived ionization corrections and abundances 
for the integrated H I and ionic column densities of the 
LLS, we would have derived a lower ionization potential, 
a lower total Ah value (« 0.3 dex), and a correspond- 
ingly higher metallicity If metal-bearing clouds consis- 
tently exhibit lower metallicity in higher Ahi components 
(Figure 7), then one may systematically overestimate the 
metallicities of the LLS by analyzing integrated column 
densities. It is also possible that such a systematic ef- 
fect afflicts current studies of the predominantly ionized 
SLLSs (Pcroux et al. 2007; Prochaska et al. 2006). In a 
predominantly ionized gas, it is of primary importance 
to precisely constrain the ionization state of the compo- 
nent that dominates the H I column density (here sub- 
system C) because this gas may dominate the total Ah 
value. 

Another important result presented here relates to us- 
ing the 0°/H° ratio to assess the metallicity of LLSs. 
It has been conventional wisdom in the ISM and QAL 
communities that the 0°/H° ratio is an excellent pre- 
dictor of O/H abundance because: (i) charge exchange 
reactions between O and H° imply similar ionization 
fractions for these atoms (i.e. O /O « H°/H); and (ii) 
oxygen contributes significantly to the metal mass den- 
sity of chemically enriched gas. The first point is valid 
even for gas subjected to an intense ionizing radiation 
field but only for radiation with a soft spectrum (stellar 
dominated). A hard radiation field (e.g. quasar domi- 
nated) is predicted to overionize O with respect to H° 
because O has a higher cross-section to hard photons 
(Sofia & Jenkins 1998). As Figure 6 indicates, a highly 
ionized gas subjected to a hard radiation field (e.g. an 

14 We expect that the inclusion of other H I clouds related to 
the LLS, e.g. the metal-poor components in subsystem B, would 
not significantly modify this result for they will have significantly 
lower JVh value. 



LLS subjected to the EUVB) should exhibit an 0°/H° 
ratio that underpredicts the true O/H value. This be- 
havior is evidenced in subsystem C of the LLS where we 
report the positive detection of Si II 1304 while no ab- 
sorption from O I 1302 is measured. In this case, the 
0°/H° ratio may undcrprcdict the O/H abundance. In 
subsystem B, however, we tentatively report the positive 
detection of O I 1302 which suggests [O/Si] > +ldex if 
we assume a quasar +galaxy radiation field. This large 
O/Si ratio suggests the gas is photoionized by a softer 
ionizing spectrum, which would require an intense, local 
starburst. We note again, however, that we cannot iden- 
tify any such source in our images of this field (O'Meara 
et al. 2006). In future papers, we will explore further the 
validity of the 0°/H° ratio as a metallicity estimator. 

Before concluding, we place the metallicities of the 
LLSs toward PKS2000-330 within the context of the 
cosmic metal budget of the z ~ 3 universe. The Ah- 
weighted mean metallicities of the LLS at z ~ 3.55 and 
the SLLSs at z w 3.19 are < [Si/H] >= -1.66 ±0.25 and 
< [Si/H] >= -1.64 ± 0.30 respectively. These central 
values are in good agreement with the mean metallicity 
of DLAs at z ss 3 (Prochaska et al. 2003), i.e. the neu- 
tral, atomic ISM gas of high z galaxies. The total Ah 
of our systems, meanwhile, is comparable to the median 
Ahi value of the DLAs (Prochaska et al. 2005). This im- 
plies that each of the LLS contributes as many metals as 
a typical DLA system. The key difference, of course, is 
that the LLSs are roughly 12 x more common than DLAs 
(Paper I). If the LLSs studied here are characteristic of 
the full LLS population (a wild speculation), the LLSs 
would contain ss 12 x the mass in metals as the DLAs. In 
this case, the LLSs would represent the dominant reser- 
voir of metals in the universe. Taking this speculation to 
the limit, we conclude that the majority of metals in our 
universe lie at the galaxy/IGM interface. This is anal- 
ogous to the notion that the metals in galaxy clusters 
resides primarily in the intracluster medium. 

6. SUMMARY 

In this manuscript, we have analyzed the Ahi value, 
ionization state, and chemical abundances of the subsys- 
tems comprising a Lyman limit system (LLS) at z ss 3.55 
toward PKS2000-330 (Table 5). We demonstrate that 
the Ahi value of an absorber whose Lya, Ly/3, and Ly7 
lines lie on the flat portion of the curve-of-growth can be 
constrained to within a few tenths dex provided analy- 
sis of the Lyman limit and the full Lyman series when 
resolved by echelle observations. As important, one can 
resolve the LLS into subsystems with distinct kinemat- 
ics and perform photoionization and abundance analyses 
for gas on (presumable) scales of several to tens kpc. 
The subsystems exhibit a presumed significant metallic- 
ity variation (> ldex) indicating incomplete mixing of 
metals in the gas comprising a galactic-scale structure. 
Finally, the following Appendix provides measurements 
of the physical properties for the SLLSs at z ss 3.19 along 
the quasar sightline. Future papers will present the Ahi 
and abundance measurements for a large sample of z ~ 3 
LLSs. 

This paper includes data gathered with the 6.5 meter 
Magellan Telescopes located at Las Campanas Obser- 
vatory, Chile. GEP and JXP were supported by NSF 
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Fig. Al. — Modeled results for the abundances (solid curves) and physical conditions (non-solid) for the gas in the cooler phase of a 
2-phase medium; see the text for details. For this cool phase, we assume photoionization equilibrium and that this phase contributes the 
fraction of C IV column density indicated on the x-axis. The results for log(/Qj^) = correspond to the single-phase scenario described 
earlier in the text (§ 4). For the scenarios where the cool phase contributes less than 1% of the observed C IV, the density of the cool 
phase must exceed njj > 0.01 cm" 3 giving clump sizes less than lOpc. For Subsystem C, we would predict anomalously low O/C ratios. 

grant AST-0307408. JO and SB acknowledge support from NSF grant AST-0307705. 

APPENDIX 

A TWO-PHASE SCENARIO FOR THE LLS AT Z w 3.55 

The results presented in the paper on the LLS atzw 3.55 were derived using a simplistic, single-phase photoionization 
model. This model assumes all of the gas is co-spatial, has constant density, has similar temperature, and experiences 
the same radiation field. The data, however, offer some indications for a more complex, multi-phase medium. First, 
the profile fits presented reveal systematically larger 6-values for the high-ion gas than the low-ion gas. This difference 
may indicate that some of the high-ion gas occurs in a separate phase from the low-ion gas, even though they have the 
same central velocity. We caution, however, that we can construct a model where the high-ion gas which coincides in 
velocity with the low-ion gas has identical 6-value. The other observation that suggests a multi-phase medium is the 
positive detection of O I absorption in subsystem B and the correspondingly large observed 0°/Si + ratio. As described 
in the text, this is difficult to accommodate within a single-phase model if we assume the intrinsic O/Si ratio has a 
roughly solar relative abundance. 

With these motivations in mind, we performed the following calculation for subsystems B and C. We adopted a 
two-phase model (cool and warm) where the cool phase is presumed to be a photoionized gas that gives rise to all 
of the observed low-ions but only a fraction /§iv °f t ne high-ion gas. The warm phase contributes the remainder 
of the high- ions observed and is otherwise ignored in this analysis. In Figure Al, we summarize the implications of 
this two-phase model as a function of /civ- The results at log/£jy = correspond, by definition, to the single-phase 
models considered in § 4. 

The solid curves in Figure Al show the chemical abundances of O, C, and Si in the cool phase. We note a strong 
dependence of O/H on /£jv which emphasizes the fact that the observed 0°/H° ratio is not a robust measure of O/H 
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in highly ionized systems. In contrast, the Si/H and C/H abundances (which are derived from the Si + /H° and C + /H° 
ratios) have only mild dependence on /civ- At the lowest /civ values considered, which correspond to \ogU w —4 in 
the cool phase, we observe for subsystem B that the O, Si, and C abundances have roughly solar relative and absolute 
abundances. The comparison cannot be well made for subsystem C as the O abundances are strictly upper limits due 
to the non-detection of O I 1302. We do note, however, that very low values for /civ imply super-solar C/O ratios. 

The gray dotted, dashed, and dot-dashed curves in Figure Al show the physical conditions of the cool phase. We 
find that the neutral fraction (1 — x) of the gas remains low (< 10%) for the range of /ciy values considered. The 
densities plotted correspond to a lower limit on the gas density where we have converted the ionization potential of 
the gas into a density by adopting the mean intensity of the EUVB atzsi 3.55 (our input spectrum from CUBA has 
log<I>EuvB = 5.6). If there is an additional, local radiation field then the densities will be higher than those plotted. 
At the lowest /civ values we find hh ^ 0.1 cm~ 3 which would be a surprisingly high value for LLS systems although 
we have little to compare against. More importantly, the inferred sizes I (plotted in units of 100 pc) for /g°°J < 10-2 
are of the order of 10 pc or smaller. In contrast, the single-phase solutions imply sizes of the order 1 to 10 kpc for the 
two subsystems. 

If we were to adopt a two-phase scenario for subsystems B and C where the majority of high-ions arise in the warm 
phase, then we are driven to a rather dense medium with very small physical size. The odds of intersection such a 
system would be extremely small unless there were millions (even billions) of such 'clouds' for every high z galaxy. For 
this reason (and the anomalous C/O and even Si/O abundances for subsystem C), we do not favor such models but 
these cannot be formally ruled out by the data. 

Finally, we comment on the implications of two-phase scenarios on the principle conclusions of this paper regarding 
the LLS system. Our conclusions on the ionization state (highly ionized) and the large difference in metal abundances 
between subsystems B and C remain unchanged. Furthermore, we would still predict a mean metallicity for the cool 
gas in the LLS of s=s 1/30 solar abundance. In short, the only major difference between the two-phase and single-phase 
scenario is in the inferred size of the absorbing gas. 

ANALYSIS OF THE SUPER LYMAN LIMIT SYSTEMS AT Z w 3.19 

Ionic Column Densities 

Strong Lya absorption at A w 5080 A is evident in the spectrum of PKS2000— 330 (Figure 1). This absorption was 
identified and analyzed in Paper I where we modeled the absorption as three Lya profiles corresponding to three 
distinct metal-line complexes at z 3.19. Each of these Lya profiles has H I column densities that likely satisfy the 
A H i criterion of a super Lyman limit system (SLLS; A H i > 10 19 cm~ 2 ). The full SLLS criteria introduced in Paper I, 
however, includes criteria one that one group all Lya lines lying within 300 km s _1 of one another into a single SLLS. 
We introduced this criterion because fits to Lya absorption lines with velocity separation Sv < 300 km s" 1 can be 
highly degenerate. In the following, however, we will treat the three fitted Lya lines as separate SLLSs and refer to 
them as SLLS a , SLLSf,, and SLLS C . We have done this primarily to investigate metallicity variations among individual 
Lyman limit systems but also because the degeneracy between the two closest Lya lines (b,c) is not too severe. We 
caution the reader, however, that SLLSf, and SLLS C should not be considered as separate in any statistical analysis of 
SLLS, especially one that draws upon the results of Paper I. 

Figure Bl presents a subset of the metal-line transitions for the three SLLSs. We note progressively stronger low-ion 
absorption from SLLS a to SLLS C . We have measured ionic column densities for these transitions using the AODM 
and present the results in Table Bl. Following Paper I (and standard practice), we fit the Lya absorption with 
three Lya lines centered at the peak optical depth of each metal-line complex: z a — 3.17292, Zb = 3.18779, z c = 
3.19169. Each line exhibits strong damping wings (Figure B2), which provides a relatively precise constraint on the 
Am values (Table Bl) independently of the assumed Doppler parameters. Our Am values are fully consistent with 
those reported in Paper I. Regarding the Doppler parameters, we have only loose constraints from the Ly/3 profiles: 
b a < 35 km s _1 , bb < 45 km s _1 , b c < 30 km s _1 . The adopted solutions are overplotted in Figure B2. As noted above, 
there is some degeneracy between the Am values of SLLSf, and SLLS C . The Lya and Ly/3 data together can be used 
to place an upper limit on the Am column density of SLLS C by assuming there is no significant Lya absorption by 
SLLSf,: A ffl < 10 193 cm~ 2 . Adopting this constraint, we set a lower limit AT m > 10 192 cm~ 2 to describe the Lya 
absorption at A > 5100A. Furthermore, we derive the constraint log A^j > log A^j — 0.1. 

Photoionization Modeling 

We have performed a similar analysis for the SLLSs at z « 3.19 toward PKS2000— 330, adopting Cloudy models with 
larger Am values and a EUVB radiation field at the appropriate redshift. Figure B3 summarizes the analysis. Note 
that for SLLS we have adopted upper limits to the C +3 and Si +3 column densities by integrating the apparent optical 
depth over the same velocity intervals as the observed low-ion absorption. There is additional high-ion absorption 
at Sv w —1225 km s _1 relative to z = 3.190 (Figure Bl) which must be arising in a more highly ionized cloud. We 
assume this gas does not contribute significantly to the observed H I absorption. In contrast to these high-ions, we do 
detect strong Si III 1206 absorption which traces the line-profiles of the low-ion transitions. Attributing this gas to 
the same phase as the low-ions, we find \ogU a = —2.8 ± 0.2 dex. If we have overestimated the ionization parameter 
by making this association, then we will have underestimated the [Si/H] abundance by w 0.3 dex. The [O/H] estimate 
(based on 0°/H°), however, establishes that SLLS has a metallicity below 1/100 solar. 
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Fig. Bl. — Metallic absorption for the SLLSs at z 3.19, referred to in the text as SLLS a , SLLSj,, and SLLS C (from left to right). All 
velocities in the figure are relative to an arbitrary z = 3.190. Blends with coincident absorption (e.g. Lyo lines) are denoted as (orange) 
dashed spectra. 

The absence of significant high-ion absorption at any velocity near SLLSb implies log Ub < —4. To be conservative, 
we adopt this limit as the central value in the following analysis noting that lower values give very similar results. 
Finally, we estimate logt/ c = —2.9 ± 0.2 dex based on the observed ionic ratios of C, Si, and Al. The absolute and 
relative abundances of the SLLSs are reported in Table B2. 
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TABLE Bl 

IONIC COLUMN DENSITIES FOR THE SLLS AT z = 3.19 



Ion 


Arest 


log/ 




log JVaodm log Af adopt 




(A) 




(km s- 1 ) 




SLLSa 



HI 


1215 


6701 


-0.3805 








19.75 ±0.15 


CII 


1334 


5323 


-0.8935 


[-1289, 


-1234 


13.88 ± 0.03 


13.88 ±0.03 


CIV 


1548 


1950 


-0.7194 


[-1289, 


-1229 


< 12.81 


< 12.81 




1550 


7700 


-1.0213 


[-1289, 


-1229 


< 13.14 




OI 


1302 


1685 


-1.3110 


[-1279, 


-1239 


14.15 ± 0.03 


14.15 ±0.03 


Al II 


1670 


7874 


0.2742 


[-1289, 


-1229 


11.59 ± 0.10 


11.59 ±0.10 


Al III 


1854 


7164 


-0.2684 


[-1289, 


-1229 


< 11.90 


< 11.90 




1862 


7895 


-0.5719 


[-1289, 


-1229 


< 12.24 




Si II 


1190 


4158 


-0.6017 


[-1289, 


-1229 


12.84 ±0.12 


12.89 ±0.05 




1193 


2897 


-0.3018 


[-1289, 


-1229 


12.91 ± 0.05 






1304 


3702 


-1.0269 


[-1289, 


-1229 


< 13.14 






1526 


7066 


-0.8962 


[-1289, 


-1229 


< 12.65 




Si III 


1206 


5000 


0.2201 


[-1299, 


-1209 


13.16 ± 0.03 


13.15 ±0.03 


Si IV 


1402 


7700 


-0.5817 


[-1289, 


-1229 


< 12.73 


< 12.73 


Fell 


1608 


4511 


-1.2366 


[-1279, 


-1229 


< 12.86 


< 12.86 



SLLSj, 



HI 


1215 


6701 


-0.3805 






19.80 ±0.15 


CII 


1334 


5323 


-0.8935 


-231, -101 


> 14.62 


> 14.62 


CIV 


1548 


1950 


-0.7194 


-251, -101 


< 12.59 


< 12.59 




1550 


7700 


-1.0213 


-251, -101 


< 12.88 




N V 


1242 


8040 


-1.1066 


[-211, -81] 


< 13.64 


< 13.64 


OI 


1302 


1685 


-1.3110 


-251, -201 


> 14.26 


> 14.26 


Al II 


1670 


7874 


0.2742 


-251, -101 


12.51 ± 0.03 


12.51 ±0.03 


Al III 


1854 


7164 


-0.2684 


-251, -101 


< 12.13 


< 12.13 




1862 


7895 


-0.5719 


-251, -101 


< 12.39 




Si II 


1193 


2897 


-0.3018 


-251, -101 


> 13.83 


13.97 ±0.03 




1526 


7066 


-0.8962 


-251, -101 


13.97 ±0.03 






1808 


0130 


-2.6603 


-251, -101 


< 14.47 




Si III 


1206 


5000 


0.2201 


-251, -101 


< 13.00 


< 13.00 


Si IV 


1393 


7550 


-0.2774 


-251, -101 


< 12.03 


< 12.03 




1402 


7700 


-0.5817 


-251, -101 


< 12.41 




Fell 


1608 


4511 


-1.2366 


-251, -101 


13.69 ± 0.04 


13.69 ±0.04 


SLLS C 


HI 


1215 


6701 


-0.3805 






19.10 ±0.15 


CII 


1334 


5323 


-0.8935 


[-18,222] 


> 14.90 


> 14.90 


CIV 


1548 


1950 


-0.7194 


[-18,222] 


13.92 ±0.03 


13.92 ±0.03 




1550 


7700 


-1.0213 


[-18,222] 


13.92 ±0.03 




N V 


1238 


8210 


-0.8041 


[22, 142] 


13.27 ±0.06 


13.27 ±0.06 


Al II 


1670 


7874 


0.2742 


[-48, 222] 


12.88 ± 0.03 


12.88 ±0.03 


Al III 


1854 


7164 


-0.2684 


[-18,222] 


12.48 ± 0.08 


12.48 ± 0.08 




1862 


7895 


-0.5719 


[-18,222] 


< 12.49 




Si II 


1190 


4158 


-0.6017 


[-18,222] 


14.14 ±0.03 


14.13 ±0.03 




1193 


2897 


-0.3018 


[-18,222] 


> 14.15 






1304 


3702 


-1.0269 


[-18,262] 


< 14.36 






1526 


7066 


-0.8962 


[-18,222] 


14.13 ±0.03 




Si IV 


1393 


7550 


-0.2774 


[-18,222] 


13.30 ± 0.03 


13.31 ±0.03 




1402 


7700 


-0.5817 


[-18,222] 


13.34 ± 0.03 




Fell 


1608 


4511 


-1.2366 


[-18,222] 


13.49 ± 0.07 


13.49 ±0.07 



a Velocity interval for the AODM relative to z = 3.190000. 



Vogt, S. S., et al. 1994, in Proc. SPIE Instrumentation in Wolfe, A. M., & Prochaska, J. X. 2000, ApJ, 545, 591 

Astronomy VIII, David L. Crawford; Eric R. Crainc; Eds., 

Volume 2198, p. 362, 362-± 
Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARA&A, 43, 
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Fig. B3. — These curves show the predictions from a series of Cloudy photoionization models for the column density ratios Si, C and 
Al ions. The models are parameterized by the ionization parameter U which corresponds to the ratio of ionizing photons to hydrogen 
atoms per unit volume. In each panel, we have tuned the Cloudy models to match the observed iVui values and metallicities of SLLS a , 
SLLSf,, and SLLS C at z £b 3.19. Overplotted on the curves are the observational constraints. We have used these to estimate the ionization 
parameter for each subsystem (designated by the shaded regions in each panel). 
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TABLE B2 
ELEMENTAL ABUNDANCES 
FOR SLLS AT z k, 3.19 



Ion [X/H] [X/Si+] 



SLLSa 



c+ 


-2.47 ±0.10 


+0.16 ±0.10 


o° 


-2.26 ±0.10 


+0.37 + 0.10 


A1+ 


-3.04 ±0.12 


-0.40 + 0.10 


Si+ 


-2.64 ±0.10 




Fc+ 


< -2.43 


< +0.20 


SLLS 6 


C+ 


> -1.64 


> -0.22 


o° 


> -2.20 


> -0.79 


A1+ 


-1.72 ±0.10 


-0.30 ±0.10 


Si+ 


-1.42 ±0.10 




Fc+ 


-1.61 ±0.10 


-0.19 ±0.10 


SLLS C 


C+ 


> -1.21 


> +0.02 


A1+ 


-1.58 ±0.11 


-0.35 + 0.10 


Si+ 


-1.24 ±0.10 




Fc+ 


-1.42 ±0.10 


-0.19 ±0.10 



Note. — In all cases, we have 
assumed a minimum error of 0.1 dcx 
due to systematic errors in the pho- 
toionization modeling. 

a Assumes a Cloudy photoionization 
model with N m = 10 19 8 cm -2 , 
[M/H] = -2.00 and log U = -2.8 ± 
0.2 dex 

b Assumes a Cloudy photoionization 
model with N m = 10 19 ' 8 cm" 2 , 
[M/H] = -1.50 and log U = -4.0 ± 
0.2 dex 

c Assumcs a Cloudy photoionization 
model with N m = 10 19 cm" 2 , 
[M/H] = -0.50 and log U = -2.9 ± 
0.2 dcx 



TABLE B3 

SUMMARY OF PROPERTIES FOR THE SLLS AT z ~ 3.19 



Property a b c 

log (N m /cm- 2 ) 19.75±0.15 19.80 ±0.15 19.10 ±0.15 

log U -2.8 ±0.2 -4.0 ±0.2 -2.9 ± 0.2 

tag(l-*)° -0.66t°-g -0.09+° °02 -1.22t° f 9 

log(Nn/cm- 2 ) 20.41 + 0.23 19.89 ±0.15 20.32 ± 0.25 

[Si/H] -2.64 ±0.16 -1.42 ±0.15 -1.20 ±0.17 
[O/H] -2.26 ± 0.25 > -2.20 

[C/H] -2.47 ±0.15 > -1.64 > -1.24 

[Fe/H] -1.61 ±0.15 -1.40 ±0.17 

Note. — Abundances assume the photoionization models as 
described in Table B2. 
a The ionization fraction x is defined as H+/H. 



